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Recently, Cormier and Holman pointed out that fluctuations of an inflaton field (f> are significantly 
enhanced in the model of spinodal inflation with a potential V(cj>) for which the second derivative 
V (2) ((f>) chan ges sign. As an application of this model, we investigate particle production in the 
, natural inflation model with a potential V(cj>) — m [1 + cos(0//)] by making use of the Hartree 

■ approximation. For typical mass scales / ~ m p i ~ 10 19 GeV, and m ~ mGUT ~ 10 16 GeV, we 

find that growth of fluctuations relevantly occurs for the initial value of inflaton 0(0) < 0.1m p i. 
Especially for 0(0) < 10 -6 m p i, maximum fluctuations are so large that secondary inflation takes 
$_j ' place by produced fluctuations. In this case, the achieved number of e-folding becomes much larger 

than in the case where an effect of spinodal instability is neglected. 
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! I. INTRODUCTION 

> 
ON 

The concept of inflationary cosmology is very attractive to describe the early stage of the universe (!]]. It not 



only solves horizon, flatness, and monopole problems the standard big bang cosmology has to face, but also provides 
large-scale density perturbations required for the structure formation 0. Inflation takes place while a scalar field 
4> called inflaton slowly rolls down toward the minimum of its potential V(<p). The inflationary period ends when 
the inflaton field begins to oscillate around the minimum of its potential. The elementary particles can be efficiently 
produced by a nonperturbative process called preheating pfl in the oscillating stage of inflaton. Then these particles 



X 

decay to other lighter particles and thermalize the universe. 



So far, there are various models of inflation. These different kinds of models can be classified in the following 
way || . The first class is the "large field" model, in which the initial value of inflaton is large and rolls down toward 
the potential minimum. Chaotic inflation || is one of the representative models of this class. The second class is the 
"small field" model, in which inflaton is small initially and slowly evolves toward the potential minimum at larger 
values of (j>. New inflation R] and natural inflation || are good examples of this model. In the first model, the 
second derivative of the potential V^ 2 \4>) usually takes positive values, but in the second model, V^{4>) can change 
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sign during inflation. The third one is the hybrid inflation model J9|, where the inflaton field has a large amount of 
potential energy at the minimum of its potential, while the vacuum energy is almost zero at the end of inflation in 
the first and second inflation models. 

Recently, Cormier and Holman JlCj ] have pointed out that fluctuations of inflaton can grow nonperturbatively during 
an inflationary stage in the second model of inflation when V^ 2 \<j)) is negative. This idea is remarkable in the sense 
that particles are effectively produced by negative instability even in the slow rolling stage of inflaton. They call this 
kind of model spinodal inflation, and investigated the nonperturbative evolution of the inflaton field making use of 

the Hartree approximation in the toy model with a potential V(cj>) = %r ^m 2 </> 2 + jr^ 4 - It was found that low 

momentum modes of fluctuations are mainly enhanced, and the evolution of the system can be described effectively 
by two classical scalar fields. Especially when <f> is close to zero initially, the amount of the maximum fluctuation 
becomes so large that secondary inflation occurs by the effect of produced particles. It was suggested that ordinary 
prediction of scale invariance of density perturbations generated during inflation would be modified by taking into 
account this spinodal effect. 

As for one example of spinodal inflation, we draw attention to the natural inflation model which was originally 
proposed by Freese et al. ||. This model is characterized by pseudo Nambu-Goldstone bosons (PNGB) which appear 
when an approximate global symmetry is spontaneously broken. The PNGB potential is expressed as V((f>) — 
m [1 + cos(</>//)], where / and m are two mass scales which characterize the shape of the potential. Considering 
PNGB as the candidate of inflaton, / and m are constrained by the requirement of sufficient inflation and primordial 
density perturbations observed by the Cosmic Background Explorer (COBE) satellite. In the case where the effect of 
spinodal instability is neglected, these mass scales are found to be / ~ m p \ ~ 10 19 GeV, and m ~ togut ~ 10 16 GeV 
respectively. While other inflation models require extremely weak coupling A in order to satisfy the constraint of 
density perturbations (for example, in the chaotic inflation model with self coupling potential, A < 10~ 13 ), the PNGB 
inflation model is preferable in the sense that two mass scales arise naturally in particle physics models. Furthermore, 
this model has an advantage in the analysis of the spinodal effect. When fluctuations of inflaton grow significantly, 
higher order terms of fluctuations play a relevant role for the evolution of the system. As compared with other more 
complicated spinodal models such as new inflation, we can handle these higher order terms in an analytic way in 
the natural inflation model. What we are concerned with is how the efficient particle production during the natural 
inflation would modify the dynamics of the system. As we will show later, secondary inflation due to fluctuations 
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pointed out in Ref. [|T0|] appears and the evolution is drastically changed if the initial value of inflaton is close to zero. 

This paper is organized as follows. In the next section, basic equations based on the Hartree approximation are 
introduced in the natural inflation model. In Sec. Ill, we study how the fluctuation of inflaton is generated during 
the inflationary stage. We will show that this effect can significantly alter the dynamics of inflation. We present our 
discussions and conclusions in the final section. 

II. BASIC EQUATIONS 



The model we consider is 



C = 



-fJ 



(2.1) 



where k 2 /8tt = G = m p 2 is Newton's gravitational constant, R is a scalar curvature, and is a minimally coupled 
inflaton field. In this paper, we adopt a potential which is the so called natural inflation type 



V(4>) = m A 



I + cos 



(2.2) 



where two mass scales m and / characterize the height and width of the potential, respectively. The typical mass 
scales are of the order / ~ m p i ~ I0 19 GeV and m ~ toqut ~ 10 16 GeV for the success of the scenario of the natural 
inflation Q. The inflaton field is initially located in the region of < 0(0) < irf, and inflation takes place while 
inflaton evolves slowly toward the minimum of its potential at <f> = nf. In order to obtain the sufficient inflation 
by which the number of e-folding exceeds N ~ 60, the initial value of inflaton is required to be close to <f> = 0. For 
example, in the case of / = I0 19 GeV, we need 0(0) < 0.5m p i. The value of <j> when inflation ends (= (j>p) depends 
on the scale /, but 4>f is close to the value irf for the typical values of / ~ m p \ ||. When inflaton begins to oscillate 
around <j) = nf, the system enters a reheating stage. 



For the potential (2.2), we find that V^{<p) is negative when inflaton evolves in the region of < < nf/2. 
This leads to the enhancement of fluctuations by spinodal instability in the realistic initial value of <p. After the 
inflaton field passes through <fi — nf/2 where V^(<j>) changes sign, fluctuations of inflaton no longer grow. Hence 
an important point for the development of fluctuations is the initial value of inflaton. The mass scales / and m also 
affect the evolution of the system. 

Let us obtain basic equations in the natural inflation model. We adopt the flat Friedmann-Robertson- Walker metric 
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ds 2 = -dt 2 + a 2 (t)dx 2 , (2.3) 

where a(i) is the scale factor, and t is the cosmic time coordinate. 

We decompose the quantum scalar field cj>(t, x) into its expectation value 4>o(t) and the quantum fluctuation S(f>(t, x) 

as 

0(t,x)=^,(t) + 50(t,x) > (2.4) 

with 

tfo(t) = Wt,x)> = ^|, (2.5) 
where /?(t) is the density matrix which satisfies the Liouvillc equation: 

i^- = [H,p(t)]. (2-6) 
TL is the time dependent Hamiltonian. The initial condition for the density matrix should be chosen to describe a 



local thermal equilibrium state. Given the initial condition, the time evolution of p(t) is known from Eq. (2.6). The 
expectation value of the energy-momentum tensor is evaluated from the time varying density matrix. Our approach 
is to solve the semiclassical Einstein equations where the source of the gravitational field is the expectation value of 
the energy-momentum tensor . 

The system we consider is out-of-equilibrium and nonperturbative, and there are some approximations which are 
suitable to describe such a state. In this paper, we adopt the Hartree mean field approximation of the nonequilibrium 
quantum field theory. This is basically a Gaussian variational approximation to the time dependent density matrix. 
Related with this issue, several authors p]"|l3| considered the large N approximation which can deal with contributions 
beyond leading order. Although it is of interest to examine the differences between two approximations, the analysis 
based on the large N approximation is left for the future work [J14J , 

Performing the Hartree factorization 

(2.7) 



,2n+l ( 2 " + 1 ) ! /s.^n 



2 n nl 

and making use of the tadpole condition 
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(<^(t,x))=0, 

the expectation value of the potential V(4>o + can be written as 

oo 1 

(V(<t> + 64)) = W 2 >"^ (2n) (0o), 

n=0 

where V^(0) = 5 n V {4>) / 5(f) n . Then the equation of the 0o field yields 

oo 1 

0o + 3i/0 o + £ ^<«0 a > n V (an+1) (0o) = 0, 



(2.8) 



(2.9) 



n=0 



(2-10) 



where a dot denotes a derivative with respect to the cosmic time coordinate, H = a/a is the Hubble parameter. 
Expanding the 5<p field by the Fourier modes as 



we obtain the following equation for the fluctuation: 

OG 



(2-11) 



a 2 ' ^ 2 n n\ 

n=0 



2n+2) 



5<f>k + 3H5(f> k + 
where the expectation value (6(f) 2 ) is represented by 

The evolution of the scale factor is written as 



= o, 



(2.12) 



(2-13) 
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^(4 2 > + ^2<W) 2 } + E W 2 )"^( 2 ")(0o) 

~" n— 



(2.14) 



where (5</> 2 ) and ((V<50) 2 ) are expressed by 



l2) = 2^ I fc2 !^| 2dfc ' 



((V^) 2 



2tt 2 



fc 4 |50 fe | 2 dfc. 



(2.15) 



(2.16) 



The quantities of Eqs. (2. IS), ( 2.15| ), and (2.16) need to be regulated in order to remove the divergences of integrals. 



Several authors considered renormalizations by the method of adiabatic regularization J12| , [l3 15 16J] and dimensional 
regulalization Jl7| . The former scheme is based on introducing a large momentum cutoff and subtracting the leading 
adiabatic orders of the fluctuation terms. The latter is the covariant regularization in which the counter terms do not 



depend on the initial state. However, this dimensional regulalization has a shortcoming that the energy-momentum 
tensor has an initial singularity. In this paper, we make use of the scheme of adiabatic regularization as in Ref. fll2| , |l3|] , 
which is suitable for numerical computations. 

In the natural inflation potential ( |2.2| ), Eqs. (2.1C ), (2.12), and ( |2.14 ) can be rewritten as 



6 + 3£T0 O - — F((5<p 2 )) sin ( J ) = 0, 



b k + 3HS<j) k + 



^-^(W 2 >)cos(- 



= 0, 



(2.17) 



(2.18) 



T 



1 • 2 



A((V^) 2 



l + F((^ 2 ))cos ( ^ 



(2.19) 



where 



F((<5</) 2 )) = exp - 



(6£ 
2/2 



(2.20) 



Note that F(((5</) 2 )) = 1 in the case of {S(f> 2 ) = 0. When 4>q is located in the region < cj)Q < nf/2 initially, the term 



in the square bracket in Eq. ( 2.18 ) is negative for small k at the first stage of inflation. This leads to the enhancement 
of fluctuations of low momentum modes. Since 0o increases for < 4>o < nf as is found in Eq. ( |2.17|) , the term 
cos(0o//) m Eq. ( p.lS| ) gradually decreases and the growth of fluctuations terminates after (j> becomes larger than 

7T//2. 



As is found by Eq. ( 2.19 ), the evolution of inflaton can be described effectively by two homogeneous fields with the 
potential 



V(4> ,a) = m 4 



' -exp cos {J 



(2.21) 



where a = y/ (54> 2 ). This potential is depicted in Fig. 1. When fluctuations do not grow relevantly as a <C /, the <po 
field slowly rolls down toward the potential minimum <f>Q — nf almost along the 0o direction in the usual manner. The 
inflationary period ends when the <po field begins to oscillate around the minimum of its potential. At this stage, the 
potential energy becomes small and the expansion rate slows down. However, in the case where a grows to the order 
of /, the evolution of the system is drastically modified. The inflaton field moves toward the a direction rather than 
the 0o direction in Fig. 1. It reaches the flat region a > 2m p i, and secondary inflation occurs there. We investigate 
for what initial values of </>o this behavior appears in the next section. 



Before analyzing the evolution of the system, we mention the initial conditions for the fluctuation. We should 
choose a conformal adiabatic vacuum state where the density matrix represents a local thermal equilibrium, which 
means that the density matrix commutes with the initial conformal Hamiltonian. This corresponds to choose the 
mode functions 6<fik as 

$0*(Q) = J—r^ , <tyk(0) = h^fe(O) - H(0)} 54> k (0), (2.22) 

y/2Wk(P) 



with 



<M0)\ R{0) 



c 2 (0) = fc 2 + M 2 (0), M 2 (0) = -— F ( W 2 » cos (^^p j - (2.23) 

where R(0) is the initial scalar curvature, and we set o(0) = 1. 

In the present model, since u> 2 becomes negative for small k, we need to modify the initial frequencies of low 
momentum modes. Following the approach performed in Ref. p2[, we adopt the initial frequency as 



fc 2 +M 2 ( 0) 
|X 2 (0)| 



lo 2 (0) = k 2 + M 2 (0) tanh ( ,7/,/^ ) . ( 2 - 24 ) 



This initial frequency coincides with the conformal vacuum frequency ( 2.23 ) for large k, and becomes positive for 
small k. This choice of the initial frequency smoothly interpolates between large and small momentum modes. An 
alternative way is to choose the initial condition as (l3) 

2(Q) _ (k 2 + \M 2 (0)\ withfc 2 <|X 2 (0)|, , , 

LJk(U) ~ \k 2 +M 2 (0) with k 2 > \M 2 (0)| . ( b) 

Although there are some subtleties about the choice of initial frequencies, we can numerically check that these different 
choices have little effect on the evolution of the system. The qualitative properties of the system are the same in 
either case of Eq. fl2.24| ) or Eq. ( |2.25f . 



We investigate the nonperturbative evolution of {S(f> 2 } with the initial condition of Eqs. (2.22) and ( p. 24 ) as the 
semiclassical problem. 

III. PARTICLE PRODUCTION IN THE NATURAL INFLATION MODEL 

In this section, we study out-of-equilibrium dynamics due to the enhancement of fluctuations in the natural inflation 
model. Let us first consider the typical case of / = 10 19 GeV~ m p \ and m = 10 16 GeV~ 10~ 3 m p i. In order to solve 
the standard cosmological puzzles of the standard big bang cosmology, the number of e-folding 
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where tf denotes the time when the slow roll period ends, is required to be N > 60. We need initial values of inflaton 
as 0(0) < 0.5m p i to obtain N > 60 in the case where the spinodal effect is not included 0. For examples, when 
0(0) = 0.5m p i, N — 71; and when 0(0) = m p i, N = 39. One may consider that the enhancement of fluctuations would 
lead to the larger amount of e-folding and relax the constraint of cj)(0) to yield N > 60. However, this is not the case. 
Fluctuations do not grow relevantly for 0(0) > 0.1m p i. We depict the evolution of the 0o field and the fluctuation 
a for the case of 0(0) = 0.1m p i in Fig. 2. We find that the maximum fluctuation at 0o — 7r/'/2 is <7 max ~ 10 _ m v \. 
Since cr^ ax /(2/ 2 ) <C 1 and F((<50 2 )) is close to unity, the evolution of the 0o field and the scale factor are almost 
the same as in the case where the growth of the fluctuation is neglected. The 0o field evolves toward the potential 
minimum at 0o = without being affected by the back reaction effect of produced particles. Inflationary period 
ends when mt ss 4.5 x 10 4 , at which the value of 0o is 0o ~ 3.0/ B. 

In the case of 0(0) < 0.1m p i, the number of e-folding becomes larger with the decrease of 0(0) since the slow roll 
period is longer. In addition to this, fluctuations are enhanced more efficiently. In TABLE I, we show numerical 
values of N and <7 max in various cases of 0(0). The number of e-folding N' where the spinodal effect is neglected is 
also presented. We find that both N and <r max increase with the decrease of 0(0) for 10~ 6 m p i < 0(0) < 10 _1 m p i. 
The growth of fluctuations continues until the 0o field reaches 0o = 7r//2. With the decrease of 0(0), since the period 
during which the inflaton field moves in the region of (0) < becomes longer, this results in the larger amount 
of the maximum fluctuation. In the case of 10~ 6 rn p i < 0(0) < 10 _1 rn p i, numerical calculations show that er max can 
be approximately written as a function of 0(0): 

m v 1 



10- 6 TO p l 



Wpi. (3.2) 



When 0(0) > 10~ 5 m p i, cr max does not exceed 0.1m p i. In this case, since (J0 2 )/2/ 2 <C 1 in Eq. (2.20), the back reaction 
effect due to fluctuations can be neglected. Although the inflaton field moves a bit toward the a direction in Fig. 1, 
it slowly rolls down toward the minimum of its potential in the usual manner. As a result, the number of e-folding 
does not change even taking into account the spinodal effect (see TABLE I). 

On the other hand, when 0(0) < 10~ 6 m p i, the fluctuation reaches cr max > m v \. In this case, produced fluctuations 
play a relevant role for the evolution of the system. For example, let us consider the case of 0(0) = 5.0 x 10 -7 m p i. As 
is shown in Fig. 3, the fluctuation reaches the maximum value er max ~ 2.3to p i at mt = 2.5 x 10 5 , where 0o = 7r//2. 



After that, particle production terminates completely because V^ 2 \(f>) changes sign, and a decreases by the expansion 
of the universe. We can see this behavior of inflaton in Fig. 1. In this case, the inflaton field evolves toward the a 
direction rather than the </>o direction, and reaches the region around a « 2.3m p i and 0o 1.5m p i. Since this region 



of the effective potential (2.21) is flatter than the region of a ss and 0o ~ 1.5m p i, the inflaton field moves slowly for 
some time. After that, it evolves along the valley around < a /m p \ < 2.3 and </>o ~ nf, and finally arrives at the 
minimum of its potential. The amount of inflation is larger than in the case where the spinodal effect is ignored as is 
found in TABLE I, because produced fluctuations provide the additional energy density in Eq. ( 2.1S| ). 

This tendency becomes stronger with the decrease of 0(0). In Fig. 4, we depict the evolution of 0o and a fields 
in the case of 0o(O) = 3.0 x 10~ 7 m p i. In this case, the fluctuation reaches the maximum value <7 max ~ 3.8m p i. One 
of the different points from the cf>(0) = 5.0 x 10 _7 m p i case is that the inflaton field stays in flat regions for a long 
time: 1 x 10 6 < mt < 7 x 10 6 . As is seen in Fig. 1, the effective potential V(<j>o, a) in the region around a « 3.8m p i 
and <fio « 1.5ro p i is very flat. Since F((S(f> 2 } max ) is much smaller than unity, V{4>q, a) takes the almost constant value 
to 4 . The third terms in Eqs. (2.17) and ( [2.1§| ) become very small in this region (since the main contributions to 
the fluctuation are due to low momentum modes, k 2 /a 2 — > 0), and the inflaton field evolves very slowly in the flat 
region V(4>q, a) « m 4 . This results in the secondary inflation supported by fluctuations. This behavior was originally 

pointed out by Cormier and Holman in the model with a potential V(4>) = ^ ^m 2 cj) 2 + |y0 4 |u|. We can expect 

that their results based on the Hartree approximation generally hold in more complex spinodal type potentials if we 
choose initial values of inflaton close to zero. When 0(0) = 3.0 x 10 _7 to p i, the secondary inflation continues much 
longer than the first inflation driven by the potential energy around <\> — at < rat < 1.3 x 10 5 (See Fig. 5). This 
means that the number of e-folding is modified to be much larger than in the case where particle production due to 
spinodal instability is neglected. In fact, as is found in TABLE I, the number of e-folding is very large asJV = 22178. 
After the secondary inflation ends, the inflaton field rolls down toward the potential minimum around 0o = n f and 
(j = 0, after which the universe enters the reheating stage. 

With the decrease of initial values of inflaton as 0(0) < 10 _7 to p i, since the fluctuation grows as o > 4m p i, inflaton 
reaches the further flat region of the potential V(<po,a). The duration of the secondary inflation becomes very long, 
and the amount of inflation is enormously large. As long as < 0o < 7r //2, the inflaton field is to roll down toward 
the a direction, and the secondary inflation continues. However, once 4>o exceeds the value of 0o = f//2) cr begins 
to decrease toward a — 0. The secondary inflation never ends in the extreme case of 0o = as was pointed out in 
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Ref. @. 

Next, we investigate the case where to and / are changed. The mass m is constrained by density perturbations 
observed by COBE. The analytic estimation neglecting particle production due to spinodal instability ^ shows that 
to is constrained by a function of / as 



1.7xl0 16 r™ Dl f(t>(tf)\Y /2 ( 15m: 



&V2 



rXP ' ) GCV ' (3 - 3) 



where b is an overall bias factor which ranges 0.7 < b < 1.3. The term s'm((f>(tf)/2f) is typically of the order of unity. 
In the case where / is of order m p \, to ranges in the region of 10 15 GeV < to < 10 16 GeV. When / is smaller than 



m p i by one order of magnitude, to decreases significantly because the exponential term in Eq. (B.E) plays a dominant 
role to determine the mass scale. If we include the effect of spinodal instability, since this would produce the spatial 
inhomogeneity, the constraint to will be changed. In order to study this issue appropriately, we have to investigate 
the evolution of metric perturbations during inflation. In the present model, however, since metric perturbations may 
be enhanced up to the nonlinear level by spinodal instability, the first order perturbation will not give the correct 
description of physics. Although we do not consider these complex issues in this paper, it would be necessary to 
include metric perturbations for a complete study of nonperturbative dynamics. 

Let us study the growth of fluctuations by changing the scale to for the fixed value of / = 10 19 GeV~ m p \. Consider 
the case of to = 10 15 GeV^ 10 _4 m p i. When particle creation by spinodal instability is neglected, since the achieved 
number of e-folding is expressed by 



16tt/ 



2 



N = n — In 

m 2 pl [ sin(<X0)/2/) 



sm(</>(i/)/2/) 



(3.4) 



it does not depend on the scale of m. Although it takes more time to terminate inflation for smaller values of to, we 
obtain N > 60 for <fi(0) < 0.5to p i, which is the same as in the m = 10 16 GeV case. As for fluctuations of inflaton, 
since (6(f) 2 ) is normalized by the square of mass m, the achieved maximum value of a becomes smaller as to decreases 
for the same initial value of (f>. Particle production relevantly occurs for the case of <fi(0) < 10 -3 to p i. Numerical 
calculations show that the maximum fluctuation cr max in the case of 10 _8 to p i < 4>(0) < 10 _3 to p i is smaller by two 
orders of magnitude than in the case of to = 10 16 GeV (see TABLE II). Namely, we find the relation 



0(0) 



-l 



10~ 2 to p1 . (3.5) 



10- 6 to p i, 

When (f)(0) < 10~ 8 to p i, cr max exceeds the order of / ~ to p i, and the secondary inflation occurs as in the case of 
m = 10 16 GeV. The number of e-folding becomes larger than in the case where the spinodal effect is neglected, and it 
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depends on the scale of m. We found that the smaller values of 0(0) are required for the development of fluctuations 
with the decrease of m. 

Finally we comment on the case where the mass / is changed. The number of e-folding is smaller with the decrease 



of /, because the potential (2.2) becomes steeper. Consider the case of / = 5.0 x 10 18 GeV ~ 0.5m p i and m = 10 15 
GeV~ 10 _4 m p i. Even in the initial value of 0(0) = 10 _1 m p i, the number of e-folding is only N — 33. In order to 
lead to the sufficient inflation as N > 60, we require initial values as 0(0) < 10 _2 m p i. Since the inflaton field rolls 
down rapidly in the regions of V^((f>) < compared with the case of / = 10 19 GeV, the growth of fluctuations is 
slower. For example, the maximum fluctuations are cr max = 6.7 x 10~ 4 m p i for 0(0) = 10~ 5 to p i; <r max = 6.7 x 10 _2 m p i 
for 0(0) = 10" 7 m p i (see TABLE III). These values are smaller than in the case of / = 10 19 GeV and m = 10 15 GeV 
for the same initial values of 0. The fluctuation grows up to the nonlinear level for 0(0) < 5 x 10~ 9 m p i, and the 
secondary inflation also occurs in this case. For the values of / which are not much smaller than the Planck order, we 
can say that fluctuations are enhanced beyond the perturbative level and can support the total amount of inflation. 

IV. CONCLUDING REMARKS AND DISCUSSIONS 

In this paper we have investigated the evolution of an inflaton field in the presence of nonperturbative behavior 
of fluctuations due to spinodal instability in the natural inflation model. Since the second derivative V^ 2 -'(0) of the 
potential V(tf>) — to 4 [1 + cos(0//)] is negative for the values of < < 7r//2, fluctuations of inflaton can grow even 
during the inflationary phase. 

The strength of the excitation of fluctuations a strongly depends on the initial value of inflaton 0(0). For typical 
mass scales / = m p i and m = 10~ 3 to p i, we have examined the dynamics of the system in various values of 0(0) by 
making use of the Hartree approximation. For the values of 0(0) < 0.5m p i, we have sufficient inflation as N > 60 
which is required to solve cosmological puzzles of the big bang cosmology. When 0(0) < 0.1m p i, fluctuations are 
relevantly enhanced with the decrease of 0(0) because duration of spinodal instability becomes longer. Since long 
wavelength modes of fluctuations are mainly enhanced, the system can be described effectively by two homogeneous 
fields with potential ( [2.21 ) . The natural inflation model has the advantage that higher order terms of fluctuations can 



be handled in an analytic way. With the increase of a, the term F(((50 2 )) in Eq. (2.2C) decreases from unity. This 



changes the evolution of the system as is found in Eqs. (2.17)-(2.19). Numerical calculations show that the maximum 



(0(O)/1O 6 m p i) 1 m p i for the case of 10 5 rn p i < 0(0) < 10 1 m p \, which means that 
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fmax is less than ~ 0.1m p i. In this case, since cr max is smaller than the scale / by one order of magnitude, the back 
reaction effect due to particle production can be neglected. When 0(0) < 10~ 6 m p i, however, a exceeds the scale of 
/ and the dynamics of inflation are altered. Since the effective potential ( [2.21 ) is flat in the region of a > 2m p \, 



secondary inflation takes place by produced fluctuations. As compared with the case where the spinodal effect is 
neglected, the number of e-foldings become much larger. The secondary inflation continues for a long time in the case 
of <f>(0) < 10 -7 TO p i. Once inflaton exceeds the value of cj>o = 7r//2, it gradually approaches the potential minimum 
around 4>q — nf and a ~ 0, after which the inflationary period terminates. 

If we change two mass scales of m and /, we obtain the smaller maximum fluctuation with the decrease of m and 
/ for the fixed initial values of (f>(0). However, if we choose smaller values of (f)(0) which are close to zero, we find 
that fluctuations can grow beyond the perturbative level to lead to the secondary inflation. The number of e-folding 
depends on the scale of m for the fixed values of / and (f)(0) if the spinodal effect is taken into account. 

We should comment on some points. The influence due to resonant particle production on cosmic background 



anisotropies was studied by several authors in the context of preheating |18| and fermion production during inflation 
p9[ in the chaotic inflation model. In the new inflation scenario, it was found that metric perturbations of super- 
horizon modes are enhanced by spinodal instability when inflaton is initially located around (f> — [ fl3| . It is interesting 
to study how scale invariance of the Harrison-Zel'dovich spectrum would be modified in the present model by making 
use of the gauge invariant formalisms of metric perturbations 20 21]]. 



As related to the treatment of the back reaction by produced particles, we relied on the Hartree approximation, 
which is essentially the mean field approximation. In the present model, this has the advantage of being able to deal 
with higher order contributions of fluctuations analytically beyond perturbation theory. On the other hand, there 
are other approaches related to back reaction issues. One of them is 2PI formalism by Cornwall et al |p2[ . Another 
approach is to add the stochastic noise term due to quantum fluctuations to the field equation, which is based on the 
closed time path formalism J23|. These formalisms may alter quantitative details obtained in this paper especially 
when fluctuations arc enhanced significantly. 

Although we studied the natural inflation model as one example of spinodal inflation, the nonperturbative evolution 
of fluctuations which leads to secondary inflation would be expected to occur in other spinodal models. As well as in the 
"small field" models such as the natural inflation and new inflation, in which inflaton is initially small, potentials with 
spinodal instability appear in the "large field" models such as the higher curvature gravity p4| ] and the nonminimally 

12 



coupled scalar field |25| by performing conformal transformations to the Einstein frame. It is of interest how the 
dynamics of inflation are modified in these models by taking into account the effect of spinodal instability. These 
issues are in consideration. 
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Note added. 

Very recently, Cormier and Holman |2lf | considered the dynamics of the spinodal instability in the same model as 
ours. Their results are consistent with our results obtained in this paper. 
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TABLE I. The number of e-folding TV and the maximum value of the fluctuation (J max for various initial values of 4>o in the 
case of / = 10 19 GeV~ m p i and m = 10 16 GeV~ 10 _3 m p i. We also attach the number of e-folding TV' where the spinodal 
effect is neglected. TV and <7 max both increase with the decrease of (f)(0) in the case of 1(T 6 < 0(O)/m p i < 1CT 1 . When 
4>(0)/m p \ < 10 -6 , the fluctuation grows significantly as o- max /m p i > 1. This leads to secondary inflation, and the amount of 
inflation drastically increases in the case of 4>(0)/m p i < 3 x 1CF 7 . 



<X0)/m p i 
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CTmax/m p l 


1 


39 


39 
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1.1 x 10" 1 


1 x 10~ 6 


776 


734 


1.1 


5 x icr 7 


1151 


771 


2.3 


3 x 10~ 7 


22178 




3.8 



TABLE II. The number of e-folding TV and the maximum value of the fluctuation a max for various initial values of (j>o in 
the case of / = 10 19 GeV~ m p i and m = 10 15 GeV~ 10~ 4 m p i. We also attach the number of e-folding TV' where the spinodal 
effect is neglected. Compared with the case of / = m p \ and m = 10~ 3 m p i, er max is smaller by two orders of magnitude for 
the same initial values of <f> while N is the same value for 0(O)/m p i > 10~ 5 . In this case, secondary inflation takes place for 
0(0) /m p i < 5 x 10~ 9 , which is smaller than in the case of / = 10 19 GeV and m = 10 16 GeV. 
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TABLE III. The number of e-folding N and the maximum value of the fluctuation a max for various initial values of <f>a in 
the case of / = 5.0 x I0 18 GeV~ 0.5m p i and m = 10 15 GeV~ 10~ 4 m p i. We also attach the number of e-folding N' where the 
spinodal effect is neglected. Both N and cr ma x take smaller values compared with the case of / = 10 19 GeV and m = 10 15 GeV 
for the same initial value of <j). Secondary inflation occurs for initial values of <j)(0)/m p \ < 5 x ICT 9 . 
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The effective two-field potential V((f>o,a) 



4 

m 



in the natural inflation 



1 + cxp (-ffr) co s y 
model. 4>o and a are normalized by the scale /. In 
the case where growth of the fluctuation of inflaton 
is neglected, this is reduced to the one-field potcn- 



. However, when the 



tial V(<f>o) = m 4 1 + cos 

fluctuation grows significantly and evolves toward 
the a direction, the system is described by two fields 
(f>o and a. 

The evolution of </>o and a fields for the initial value 
of </>(0) = O.lmpi in the case of / = 10 19 GeV^ m pl 
and m = 10 16 GeV^ 10~ 3 77i p i. Both fields are nor- 
malized by m p \. Although fluctuations grow at the 
initial stage, the maximum value cr max ~ 10~ 5 m p i 
achieved in this case is much smaller than the value 
/. The <po field evolves toward the potential mini- 
mum <po = 7r/, after which the universe enters the 
reheating stage. In this case, the enhancement of 
fluctuations hardly affects the evolution of the <p 
field and the scale factor. 

The evolution of (f>o and a fields for the initial value 
of (f>(0) = 5.0 x lCT 7 m pl in the case of / = 10 19 GeV~ 
m p i and m = 10 16 GeV^ 10~ 3 m p i. The fluctuation 
a reaches the maximum value cr max = 2.3m p i when 
4>v = 7r//2. After that, a decreases because spinodal 
instability is absent. The inflaton field finally rolls 
down toward the potential minimum with <po = nf 
and (7 rj 0. 



The evolution of (f>o and a fields for the initial value 
of (f)(0) = 3.0 x 10" 7 m pl in the case of / = 10 19 
GeV~ m p] and m = 10 16 GcV- 10- 3 m pl . The fluc- 
tuation a reaches the maximum value cr max = 3.8to p i 
when 0o = tt//2. The secondary inflation due to 
fluctuations occurs for 1 x 10 6 < rat < 7 x 10 6 . In 
this region, since the effective two-field potential is 
very flat, the inflaton field evolves very slowly. Fi- 
nally, inflaton is trapped in the potential minimum 
(f> = nf and a w at mt = 7.6 x 10 6 . 

The evolution of the Hubble parameter H for the 
initial value of (f>(0) = 3.0 x 10~ 7 77i p i in the case 
of / = 10 19 GeV~ TO p i and m = 10 16 GeV~ 
10~ 3 to p i. The first inflation occurs around 0o ~ 

for < mt < 1.3 x 10 5 , which is followed by 
the secondary inflation caused by fluctuations for 

1 x 10 6 < mt < 7 x 10 6 . Since the duration of this 
secondary inflation is long, the achieved number of 
e-folding is very large as N = 22178. 



17 



Fig. 1 




Fig. 5 



0.005 



0.004 



0.003 



0.002 



0.001 



-0.001 



2 10 b 



4 10 6 



_1 I I L_ 



6 10 6 



8 10 6 



1 1( 



